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Abstract 

The high density effective theory recently introduced by Hong and Hsu to describe 
ultradense relativistic fermionic matter is used to calculate the tree-level forward 
scattering amplitude between two particles at the Fermi surface. While the direct 
term correctly reproduces that of the underlying gauge theory, the exchange term has 
the wrong sign. The physical consequences are discussed in the context of Landau's 
theoretical description of the Fermi liquid. 
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1 Introduction 



When studying any strongly-interacting system, an effective description which high- 
lights the low energy degrees of freedom is highly desirable. Such a programme has 
recently been initiated by Hong and collaborators for the case of systems with a large 
number density of quarks, in which case the simplified description keeps as excitations 
about the ground state just those particle and hole states in the immediate vicinity of 
the Fermi surface. The resulting high density effective theory (HDET) pQ is derived 
from QCD in very much the same way as the heavy quark effective theory (HQET) [2], 
which has proved so successful in the description of bound states of heavy and light 
quarks. In HQET, heavy quarks of mass M moving with 4-velocity t> M are described 
by fields h v defined in terms of the original quark fields by 

h v (x) = e iMv ^P v tfj(x) with P v = \(l + i,). (1.1) 

The projector P v enables the kinetic term for the h v fields to be independent of the 
large scale M, while the phase factor ensures that even when the h v fields are restricted 
to carry some residual momentum p <C M, gluons which scatter off the heavy quark 
exchange momentum with a state carrying physical momentum P ~ O(M). HDET 
employs a similar manipulation with M replaced by the chemical potential /x, which 
for a degenerate system at zero temperature may be identified with the Fermi energy. 
It is then possible to devise a theory in which both particle and hole excitations 
appear with equal measure in the path integral [Sj. 

The HDET approach has succeeded in calculating the screening mass of electric 
gluons in quark matter, and in calculating the strength of the gap induced by a 
BCS instability to diquark pairing in the color-flavor locked (CFL) channel pQ. Most 
intriguingly, in Euclidean metric it yields at leading order in A/fi (where A is the 
low energy scale of interest) a positive definite path integral measure [Sj. This has 
enabled a Vafa-Witten style proof that the CFL phase is the true ground state of 
QCD at asymptotically high density jlj, as well as raising the possibility that HDET 
might form a starting point for non-perturbative lattice simulations of dense matter. 

There is, however, one important physical respect in which dense matter differs 
from the heavy-light systems described by HQET. In the latter case the only impor- 
tant interactions of the heavy quarks are small-angle scattering by t-channel exchange 
of gluons with light quarks. In quark matter scattering takes place between identical 
particles, implying that interactions in the it-channel, in which large momenta may be 
exchanged between two non-parallel particles at the Fermi surface, are also relevant. 
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In many-body and nuclear physics these two distinct contributions are known respec- 
tively as "direct" and "exchange" processes. While HQET and HDET both capture 
the essential physics of the direct interaction, for HDET one should also check its 
behaviour in the exchange channel. 

After reviewing the version of HDET as set out in Refs. 0] in the next section, 
I will outline in Sec. E3 an explicit calculation of the amplitude for forward scatter- 
ing between two particle states at the Fermi surface, using both gauge theory and 
HDET, and show that while HDET successfully reproduces the direct term, there is a 
mismatch for the exchange. Forward scattering is interesting for degenerate systems, 
and indeed, is known to be the only non-irrelevant interaction at the Fermi surface 
in the renormalisation group sense 0. In Sec. |U I will show that this amplitude is 
a central feature in an alternative and much older phenomenological description of 
dense matter, the Fermi liquid [B]. The consequences of an incorrect description of 
exchange interactions for the detailed relations between Fermi liquid parameters such 
as the Fermi energy, momentum and velocity, as well as for the behaviour of collective 
excitations, are then set out in Sec. El Finally I present some brief conclusions. 



We begin by reviewing the derivation of the high density effective theory (HDET), 
starting from the fermionic part of the gauge theory Lagrangian density and including 
a chemical potential /i (we neglect any bare quark mass): 



where the covariant derivative = + igA^. Rewrite the physical quark 3- 
momentum as k = /ip + p where the "residual" momentum, giving the distance 
from the Fermi surface, satisfies p || k. HDET assumes |p| -C /x, and is expressed in 
terms of fields most naturally thought of in reciprocal space as functions of p. This 
is achieved via the decomposition into "fast" ip- an d "slow" ip + degrees of freedom 



2 High Density Effective Theory 






(2.2) 



where |p| = 1 and ijj± are eigenstates of the projection operators 




(2.3) 
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The matrix a = 707, and where needed we use the representation 



a 



s ={, )■ < 2 - 4 > 

We also define a new covariant derivative D^: 

D, = d, + ige-^A^v = d, + igA^. (2.5) 

In real space the fields must be constructed using the operator p = (— «/a/V 2 )V. In 
terms of the new fields the quark Lagrangian (j2.1|) becomes 

C = $+7o(l, pTi£>^+ + $_7o[(l, -P)"*^ + M-*l>- + [$-iHD ll i/> + + h.c. 

= TP + i^D^ + + V-(nf^ + 2^70)^- - 9 [^-7±A^+ + h - c ] ( 2 - 6 ) 

where we define 7^ = (70, P7-p) M , 7^ = 7* 1- 7n\ and the form of the final term follows 
from 7|| 1 <9 A1 = $ , IxPn = ®- We nave used the identities [I] 

P T rP ± = jo(i,±i>YP± ; P±YP± = ilP±. (2.7) 

For p <C p it is possible to integrate out the fast modes to yield a low energy 
effective Lagrangian written solely in terms of At tree level this is done via the 
equation of motion 

where D\\ = (1, — p^-D^. We obtain the HDET effective Lagrangian as a derivative 
expansion, ie. in powers of A/p where A <C p denotes a momentum scale of physical 
interest : 

9 j. 4 \2„,. , ^ f P J 



Chdet = V+i7f ^+ + f^+7o(7^)V>+ + (^J • ( 2 -9) 

The leading order term of HDET should therefore be able to describe the dynamics of 
low energy excitations near the Fermi surface, in which any momentum transferred by 
scattering off a gluon should be much smaller than the Fermi momentum. At higher 
order, further terms in the HDET action such as a four-fermi contact interaction and 
a gluon screening mass are generated [T]. 

In fact, HDET has to be subtly modified in order to preserve unitarity. The ?/>_ 
excitations are nothing other than anti-particles, whose physical impact in superdense 
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matter with p <C p, is relatively unimportant. However, just as in a metal or a 
semiconductor, we must also take hole excitations in the Fermi sea into account, which 
have identical quantum numbers to anti-particles, but play a much more important 
role. In order to do this consistently we modify the definitions (|2.21 12. 5|) to read [3]: 

iP(x) = e iX 7p+(x) ; A^x) = e~ iX A^x)e iX with X = /i(x.p) (a.p) . (2.10) 

Note that ip+(p) can now satisfy either of P±(p)if} + = ip + depending on whether the 
physical state it describes is a particle or a hole, in either case with momentum p. 
It is not possible, however, to derive HDET from the original gauge theory starting 
with the projection (|2.10j) . The quark-gluon vertex arising from the leading order 
term ip + i^D^ilj + in Chdet now reads 

C qqg = -g I A»(q)e iqx 4)+{p l )e- lv ' x e- i ^' d -»\ l0 ,t>-^ (2-11) 

J X 

After integrating over x we find 

C gqg = -gA"(q) [Mp')P4p') l0 (l, +P%P + (p)<Mp)5 4 (V + pi>') -(p + w)-q 
+ Mp')P + (p') l0 (l, -p),P-(p)*P+(p)S 4 ((p' ~ MP') - (P - /4>) - q 

+ $+{p')P-(p'ho{i, -i>)vP-{p)4>+{p)5 A ({p + nv' ) - (p - ^p) - 5 

+ ^ + (p')P + (p')7o(1,+p)^+(p)^+(p)5 4 ((p' -W)-{P + W)-<1)\ (2-12) 

the four terms describing respectively particle or a hole scattering off a gluon, and 
particle-hole creation and annihilation. The effect of the phase factor in ()2.10|) has 
been to ensure that the gluon scatters off states carrying physical rather than residual 
momenta. The Dirac structure in principle is given by expressions such as 

P-(p')l P+(p) = j [1 + P'-P + «-(P + P) + H55.P x p] (2.13) 

Note, however, that since by assumption p'.p = 1 + 0(p 2 //i 2 ), p'xp = 0(p/p), at 
leading order in HDET we can consistently approximate 

P^(p'hoP±(p)~loP±(p)- (2-14) 

In this same limit P±{p') r )oP±{p) ~ 0. 

In Euclidean metric, the leading term of the HDET Lagrangian ()2.9|) is of the 
form %() + M%p + where M is both anti-hermitian and anticommutes with 75, implying 
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that the functional measure of the path integral detM is real and positive. This raises 
the attractive possibility of using the HDET action as a basis for non-perturbative 
simulations of systems at asymptotically high quark densities using lattice Monte 
Carlo methods [SJ HJ. The higher order terms in (|2.9|) describing corrections of 0(A/ /i) 
are of the form ip + Hip + where H is hermitian. In general, det(M + H) is complex, 
leading to a restoration of the notorious "sign problem" at finite density. It may prove 
possible, however, to find means to treat corrections for A//i J$ 1, much as corrections 
to QCD simulations at /i = can be calculated for /x/T ^ 1 [Jj. 

3 Forward Scattering 

Let us consider a very simple process, tree-level forward scattering of two relativistic 
particles at the Fermi surface, by comparing calculations made using both HDET and 
the full underlying quantum field theory, which for simplicity we will initially assume 
to be QED. The relevant Feynman diagrams are shown in Fig.^ In Minkowski metric 




Figure 1: Feynman diagrams for forward scattering, showing direct (left) and ex- 
change (right) contributions. 

and Feynman gauge the QED scattering amplitude for the direct process is 

A dir (k, k') = —?—M dir (k, k') (3.1) 
4£ k e k / 

with 

iM dir = ie 2 ^—u(k, s^uik, s)u(k', s')j v u{k', s'). (3.2) 

Here, e is the electron charge, an infra-red regulator mass A has been introduced for 
the photon, and u, u are plane-wave spinors normalised to uu = 2m, with m the 
electron mass. The four-momentum k^ = (e^, k) M , with the energy e k of course equal 
to jj, at the Fermi surface. 



6 



To evaluate the spin-symmetric part of the amplitude we average over s, s' to 
obtain 

A dir (k, k') = -— — ^tr[ 7M (^ + m ))tr[r(k" + m)} = - — . (3.3) 

The exchange contribution includes a relative minus sign due to Fermi statistics: 

! .[_( £k£k ,_ k .k') + 2m 2 ]. (3.4) 



2e k ev (k - k') 2 + A 2 

Note that ko = k' for two states at the Fermi surface. In the limit m, A — > we have 
£k = |k| = \i and jHj 



A% r ED (k,k') = -^(l-cos0), (3.5) 



^ B d(*,*') = (3-6) 

where 9 is the angle between the two particle momenta. 

To repeat the calculation using HDET a necessary ingredient are the plane wave 
states. Eigenstates ty± of the momentum space HDET Dirac equation 

Po*±(p) = 3.pP±(p)*±(p) (3.7) 

are given by 



ty+(p, s) = ( ) with eigenvalue po = +|p| 

X° 



ty-(p, s) = ( ^p^sj with eigenvalue po = —|p|- (3.8) 
The two-spinors 0, x are given by 1 = x 2 — (^j > 2 — X 1 = ( ^ j > whence 

^± = 2 ; # ± * ± = ; ^^ ± (p, S )^ ± (p, S )=7 M (l,±p) At . (3.9) 

Note also the important identity 

*±(p,s) = * T (-p,-s). (3.10) 
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Now let's use HDET, and in particular the qqA vertex described by the first term 
of ()2.12jl . to calculate the spin-symmetric forward scattering amplitude between two 
particles at the Fermi surface. First we examine the direct term; with the normalisa- 
tion ()3.9|) we have 



16 

s,s' 

J S H(1, +p)a**+(p> s')7o(l, +P')^+(P', «') 



16 -A 2 

s,s' 

= - I ^tr[ 7 o7p]tr[7o7 CT ](l,p) P (l,p0 CT (l-RP0=^&Dfep / )- (3.11) 

The direct term calculated in HDET thus agrees with that of the full theory l)H.5jl . 
For the exchange term, 

s,s' 

(3.12) 

with q = (|p'| + /x)p' — (|p| + /i)p (as before, for two states near the Fermi surface 
go ~ 0). Performing the average over spins and resulting trace we find 

2 2 

AhdetM) = -^QfioQaO - £oo<U(l,p) p (l,py(l - Pp') = ~{1 - cos ^)- 



(3.13) 

Finally, q 2 = — 2/i 2 (l — cos#)(l + 0(p/ [/,)), so that to leading order the result is 

,2 



A% DET (p,p') = ^(l + cos9). (3.14) 

o/i 



If the calculation is repeated with the more accurate vertex factor 1)2.13)1 . then the 



answer is 



e 2 



^hdeAp,p') = 5-2 
8/r 



1 + cosfl) 2 1 



H — sin 6(1 — cos ( 



(3.15) 



2 

As anticipated, this result agrees with ()3.14)1 up to terms of 0(p 2 //i 2 ). In either case, 
though, the HDET result, though of the same order of magnitude, actually disagrees 
even in overall sign with that of QED ([3.6)1 . 

Since QED and HDET disagree even at leading order in p/fi, it is difficult to 
see how the discrepancy between (13.6)1 and ()3.14j) can be sorted out at higher order; 
indeed, the next term (2/z) _1 ?/> + 7o(e4 1 ) 2 ip+ in the derivative expansion would correct 
A e H DET only at 0(e 4 ). In fact, to reconcile the two approaches we need to step back 
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to the mixed term in ()2.6|) before ■?/>_ is eliminated: e[ip_jfL ± ip + + h.c.]. This yields 
a vertex — iejo(®, <3 — p)m> which can contribute to particle - particle scattering once 
we realise that an anti-particle created by ip-{p) is indistinguishable from a particle 
with momentum —p. The new vertex does not contribute to direct scattering, but 
does yield three additional contributions to A ex at 0(e 2 ), eg: 



e 2 



A +- = -Y^E^+^^07o(l,p)^ + (p,s)^4-P,-s)7o(0,a-p')^ + (p'^ / 



e 2 r / 1 

~V [P^f 09 ^ + 9pi9ao) + p.p(2^ p o^o - 9 P a9oo)\ (1, p) p (l, p) c 

e 2 (l + cos#) 2 



(3.16) 



8/i 2 1 — cos# ' 
where use has been made of (I3.1UJ1 . We then find that to 0(e 2 ) 

A e * + + A e *_ + A e * + + A e *_ = A% ED (3.17) 

where ++ denotes the original HDET term (|3.14|) . We conclude that once the ip- 
degrees of freedom have been eliminated to define the effective theory solely in terms 
of "0 + , the exchange amplitude (J3.6|) cannot be faithfully reproduced at tree level. 

4 The Fermi Liquid 

Why is this important? The forward scattering amplitude plays a central role in 
another, much older phenomenological description of degenerate matter, the Fermi 
liquid jHlEl- The essential physical idea is that the dominant low-energy excitations 
in the neighbourhood of the Fermi surface are quasiparticle states having energy e^, 
width of 0(£k — /i) 2 and equilibrium distribution 

= exp(4*) + l- (4 ' 1) 
For temperature T = we expect to have the form 

e k ~V + p F {\k\-k F ) (4.2) 

where k F , f3 F = |Vk£|k|=fc F | are respectively the Fermi momentum and Fermi velocity, 
which at this level are phenomenological parameters. Leading order HDET assigns 
them the free-field values k F = [i, j3 F = 1. Fermi liquid theory derives quantitative 
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power via the equation for the variation in quasiparticle energy in response to a 
departure Sn from equilibrium: 

/d 3 k' 
j—^Fkvdnv. (4.3) 

The Fermi liquid interaction .Fk,k' is given to lowest order in perturbation theory by 

maun 

^k,k' = -A{k, k') = -.A(p, p') (4.4) 

where in the last step we have assumed that the states are so close to the Fermi 
surface that the matrix element only depends on the relative orientation of the particle 
momenta. If T can be calculated in some scheme, it is then possible to obtain 
quantitative relations between the parameters //, kp and /3p, as well as derive the 
velocities of collective excitations I3IE1IH]- 




Figure 2: Schwinger- Dyson equation for the scattering amplitude —T(p',p,q). 

It may be helpful to outline the origin of the relation (|4.4|) - to my knowledge 
there is no argument more elegant than the original one presented by Landau in one 
of the earliest papers on non-perturbative field theory jTQl HI] • Consider quasiparticle 
scattering with the kinematics sketched in Fig. EJ the full amplitude is then given by 

T(p',p,q) = r°(p',p)-i J ^-r(p',£,q)S(£)S(£'(£,q))T (q,p). (4.5) 

We wish to focus on the singular behaviour of T as the momentum transfer q — > 0, 
so separate out the part T (shown as a simple point vertex in Fig. |2J) which is 
regular in this limit. The singular behaviour results from the quasiparticle propagator 
product S(£)S(£'), and in particular from the contribution of a particle-hole pair, in 
which no net quark number flows through the diagram, for loop momenta on the 
Fermi surface. Using HDET propagators which are functions of residual rather than 
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physical momenta, but with spinor indices suppressed for simplicity, this product 
contains terms like 

S(£)S(£'(£, q))=gx i — — 1 (4.6) 

where we have written uo = —go, the ie terms are chosen so that the first factor 
descibes a particle of momentum £, the second a hole of momentum —£', and the 
kinematics enforced by eg. the (^-functions of the third or fourth terms of (|2.12|) 
are such that |1'| = — 11| + |q| cos#, where 9 is the angle between q and 1 and we 
have assumed all momenta <C \i so that the curvature of the Fermi surface can be 
neglected and 1' || 1. The factor q arises from a trace over internal degrees of freedom, 
and counts the degeneracy of each momentum state. Within the integral in (|4.5jl the 
product may be approximated by 



S(£)S(£') ~ 27r z 14 , ^ (W|l|) + G(£) (4.7) 





q| 


cos 9 


to - 




|q| 


cos 6* 



where the pole structure has resulted in the first term being supported only for loop 
momenta corresponding to physical momenta on the Fermi surface, and the second 
term G is once again regular as q — > 0. The factor in the numerator arises because 
the kinematics restricts the integral over 3-momentum to < |1| < |q| cos#. 

The rest of the argument proceeds by examining the resulting integral equation 

T(p',p,q) = T°(p',p) - ij J^r(p',£,q)G(£)r (£,p) 



+ WTs f dnr( P ',£,q)T\£, P ) |qC " Sg (4.8) 
(2Tr)' i J lu — |q| cost^ 

in two distinct kinematic regimes as q — > 0. First consider \q\/u — > 0: 

r"(p',p) = lim T{p',p,q) = rV,P) - i / -^r>(j/,QG(£)T°(l,p). (4.9) 

q^O, |q|/w-»0 J (27T) 4 

In fact, it is possible to eliminate dependence on the regular functions T and G by 
substituting fll~3J) in to yield 



T(p f ,p,q) = Vtf,p) + I dnr(p',£,q)T»(£,p) |q °° S * (4.10) 

(27r) d y u; — |q| cost' 

Now we can consider the opposite limit u/\q\ — > in (|4.10|) : 

r fc (p', P ) = lim r(p', P ,g) = r"(p',p) - / ^r fc ( P ',£)r-(£,p). (4.n) 

g-*0,w/|q|-»0 (27TJ' 1 J 
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The function T k describes forward scattering for quasiparticles located on the Fermi 
surface, for which all physical scattering processes have u = 0, and can be identified 
with —A calculated in the previous section. In order to find a physical interpretation 
of r w , we consider (j4.10J) in the vicinity of a pole of T(p',p, q) considered as a function 
of q; in this case the first term on the RHS can be neglected and the argument p' which 
plays a passive role can be suppressed. At the Fermi surface the residual momenta p 
and I can be replaced by their corresponding unit 3-vectors. With the definition 

Hp) = - t ^-T(p,q) (4.12) 
u - pq 

eqn. ()4.10j) becomes an integral equation for excitations propagating with wavevector 
q and amplitude $: 

(a; - P .q)$ ( P ) = (P-q)^ / ^(1, p)$(l). (4.13) 

This equation 1 describes collective excitations of the shape of the Fermi surface known 
as zero sound, one of the principal predictions of the Fermi liquid theory ^2], with 
T w playing the role of Fermi liquid interaction T . Eqn. 1)4.11)1 thus yields the full 
relation between A and T: to lowest order in perturbation theory we recover (j4.4j) . 



5 HDET Confronts the Fermi Liquid 

There are many situations in strong interaction models, such as the chiral limit of 
the NJL model ^3], or the color symmetric channel in QCD jS], where the direct 
contribution to A and hence T vanishes. In this situation HDET would yield a Fermi 
liquid interaction of the wrong sign, with potentially unphysical consequences. 
To see how this works, define the spin-symmetric Landau parameters 

if = (2/ + 1) I ^Pi(cos9)F ky (9) (5.1) 

where the integral ranges over all angles between k and k' taken on the Fermi surface. 
I now quote without proof the following relations for a relativistic Fermi liquid [B] : 



kp gk 



2 



Pf = --^f/f; (5.2) 
dp, 2n 



2 



in - |/f - (5-3) 



dn pgkj 



If proper account is taken of the absence of full Lorentz invariance, then all 3-momenta in the 
argument leading to 1)4.13)1 are rescaled by a factor (3f- 
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Here n = Qkp/GTt 2 is the quark number density. Generalising the exchange result 
((HUH) to QCD with N c colors and Nf flavors of massless quark (note the direct term 
vanishes in the color symmetric channel) we have 



; /f = 0: S = 2N ' N >- (5 ' 4) 

Integrating ()5.3j) we thus find to 0(g 2 ): 



a - kF -i 3 2 (^-i) 

For the HDET exchange interaction ([3. 14)1 . however, the corresponding results are 

fo = fi = -9 2 iGN2 Nf li ( 5 ' 6 ) 

implying 

P*-j{ 1 + 9 ~48N^)- 1 + 9 lUN e ir* ■ (5 ' 7) 
Thus the Fermi velocity is predicted to be superluminal in this approach. 

At 0(g 2 ) we can correct for the mismatch between ()5.5|) and (J5.7|) by modifying 
the leading order HDET Lagrangian to 

^W = ^+(7o,c 1 p7-P) M ^+ with A ft = e~ iX 'A /t e iX ' ; X' = c 2 X, (5.8) 

with the constants chosen to be 

9 11(N 2 -1) JN 2 -l) , . 

However, whilst it is feasible to repair HDET to reproduce parameters directly related 
to the fundamental quanta, there are other phenomena where the failure to deal with 
the exchange interaction correctly may be more serious. For instance, in theories 
where the only non-trivial Landau parameter is /q , eqn. (j4.13j) may be solved for the 
zero sound velocity 3q = uj/\q\ [T2] : 

8 = 3 F (l + 2exp (-*£faX) . (5.10) 



The collective excitations are thus extremely sensitive to the sign of /q 
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6 Conclusion 



The main result of this paper is that if the tp_ degrees of freedom defined by the 
projection (|2.2|) are eliminated, the resulting effective theory incorrectly describes 
exchange interactions between quasiparticle states at the Fermi surface at tree level, 
as exemplified by the explicit perturbative QED calculation of Sec. |3] 

It should be stressed that this result does not invalidate the HDET programme; 
the failure to reproduce exchange amplitudes may be compensated by introducing a 
four-fermi contact term of 0(g 2 /fi 2 ), as outlined by Schafer 14 . In his approach all 
operators consistent with the symmetries of the underlying gauge theory are included 
from the start, and coefficients systematically determined by matching at a suitably 
chosen scale A ~ \i. Within the approach adopted here in which the leading order 
terms are "derived" via the projection (|2.10|) . the correction arises from integrating 
out ■?/>_ states from the loop of Fig. EJ but with the naive power counting of 0(g A / /i 2 ) 
non-perturbatively enhanced to 0(g 2 //i 2 ). 

How serious a problem is this? It is worth recalling that currently few people be- 
lieve that quark matter is a Fermi liquid; rather it is supposed to be a superconductor 
with energy gap A ~ O(10)MeV at the Fermi surface. The infra-red singularity lead- 
ing to the key relation (j4.11j) for the Fermi liquid interaction is thus cut off, and 
the dire consequences predicted in Sec. 03 no longer hold. The relativistic degenerate 
electrons in the interior of a white dwarf are expected to form a Fermi liquid, but in 
this case a photon screening mass A ~ 0(efi) should be included in the calculation, 
which implies that the exchange amplitude (jH.fi |) is only 0(a) compared to the direct 
term. 

At the very least though, these considerations suggest that careful physical argu- 
ments should be supplied before the HDET approach is applied, particularly since one 
of the main goals of any conceivable numerical implementation would be to examine 
the parameter range of its applicability. 
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